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I. INTRODUCTION 

"The problem of turbulence" is often hailed as one 
of the last open problems of classical physics. In fact, 
there is no single "problem of turbulence" ; rather, there 
are many inter-related problems, some of which had seen 
significant progress in recent years, and some are still 
open and inviting further research. The aim of this short 
review is to explain where fundamental progress has been 
made and where, in the opinion of the present writers, 
there are opportunities for further research. 



There are many ways to set a fluid into turbulent mo- 
tion. Examples include creating a large pressure gradient 
in a channel or a pipe, pulling a grid through a fluid, mov- 
ing one or more boundaries to create a high shear and 
forcing a high thermal gradient. Customarily the vigor of 
forcing is measured by the Reynolds number Re, defined 
as Re = UL/v where L is the scale of the forcing, U is 
the characteristic velocity of the fluid at that same scale, 
and V is the kinematic viscosity. The higher the Reynolds 
number the larger is the range of scales involved in the 
turbulent motion, roughly from the scale L itself (known 
as the "outer" or "integral" scale) down to the so-called 
"viscous" scale 77 which decreases as Re~'^/^ [1]. For large 
Re a turbulent flow exhibits an erratic dependence of the 
velocity fleld on the position in the fluid and on time. 
For that reason it is universally accepted that a statisti- 
cal description of turbulence is called for, such that the 
objects of interest are almost invariably mean quantities 
(over time, space or an ensemble, depending on the ap- 
plication), fluctuations about the mean quantities, and 
correlation functions defined by these fiuctuations; pre- 
cise definitions will be given below. Thus, the crucial 
scientific questions deal typically with the universality of 
the statistical objects, universality with respect of the 
change of the fluid, or universality with respect to the 
change of forcing mechanisms. We will see that this uni- 
versality issue binds together the various aspects of tur- 
bulence to be discussed below into a common quest — the 
quest for understanding those aspects of the phenomenon 
that transcend particular examples. We will strive to un- 
derline instances when this quest has been successful and 
when doubts remain. 

The structure of this review is as follows: in Sect. 
II we discuss the statistical theory of homogeneous and 
isotropic turbulence and focus on the anomalous scal- 
ing exponents of correlation functions. For a part of the 
community this represented the important open prob- 
lem in turbulence, and indeed great progress had been 
achieved here. In Sect. Ill we address homogeneous 
but anisotropic turbulence and present recent progress 
in understanding how to extract information about the 
isotropic statistical objects, and how to characterize the 
anisotropic contributions. Section IV deals with wall- 
bounded turbulence where both isotropy and homogene- 
ity are lost (being actually the norm in practice, rather 
than the exception). We focus on the controversial issue 
of the log law versus power laws, clarifying the scaling 
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assumptions underlying each of these approaches and re- 
placing them by a universal scaling function; we show 
that this achieves an excellent modeling of channel or 
pipe flows. In Sect. V we consider turbulence with addi- 
tives (like polymers or bubbles) and review the progress 
in understanding drag reduction by such additives. Sec- 
tion VI discusses problems in thermal convection, with 
emphasis on recent work. Finally, Sect. VII provides a 
selective account of the problems that have come to the 
fore in superfluid turbulence, sometimes bearing directly 
on its classical counterpart. The article concludes with a 
summary of outlook. 

II. ANOMALOUS SCALING IN 
HOMOGENEOUS AND ISOTROPIC 
TURBULENCE 

A riddle of central interest for more than half a cen- 
tury to the theorist and the experimentalist alike con- 
cerns the numerical values of the scaling exponents that 
characterize the correlation and structure functions in 
homogeneous and isotropic turbulence. Before stating 
the problem one should note that strictly homogeneous 
and isotropic state of a turbulent flow is not achievable in 
experiments; typically the same forcing mechanism that 
creates the turbulent flow is also responsible for breaking 
homogeneity or isotropy. Nevertheless, some reasonable 
approximations have been attained. To get a closer ap- 
proximation, one has to resort to numerical simulations. 
For a long time, the Reynolds number of simulations was 
limited by numerical resolution and by storage capabili- 
ties, but this situation has improved tremendously in the 
past few years. Indeed, as an idealized state of turbulence 
which incorporates the essentials of the nonlinear transfer 
of energy among scales, homogeneous and isotropic tur- 
bulence has gained a time-honored status in the history 
of turbulence research. 

Consider then the velocity field u{r, t) which satisfies 
the Navier-Stokes equations 

^+u-Vu=-Vp + vW'^u + f , (1) 
ot 

where p is the pressure and / the (isotropic and homo- 
geneous) forcing that creates the (isotropic and homoge- 
nous) turbulent flow. Deflning by (...) an average over 
time, we realize that {u) = everywhere in this flow. 
On the other hand, correlations of u are of interest, and 
we define the so-called "unfused" nth order correlation 
function r„ as 

Tn{ri,ti,r2,t2,... rn, tn)= {u{ri,ti)u{r2,t2) . . . tx(r„, i„)) 

(2) 

When all the times ti arc the same, U = t, we get the 
equal-time correlation function i^„(ri, r2, . . . r„) which, 
for a forcing that is stationary in time, is a time- 
independent function of the n(n — 1)/2 distances between 
the points of measurements, due to homogeneity. An 



even more contracted object is the so called "longitudi- 
nal structure function" 5'„, 

Sn{R) = {{[u{r + R,t)- u{r, t)] ■ R/RD , (3) 

which can be obtained by sums and differences of corre- 
lation functions i^„, together with some fusion of coor- 
dinates [2]. On the basis of evidence from experiments 
and simulations, it has been stipulated (although never 
proven) that Sn is a homogeneous function of its argu- 
ments when the distance R is within the so-called "iner- 
tial range" t] <^ R<^ L in the sense that 

Sni\R) = X^"Sn{R) . (4) 

The central question concerns the numerical values of 
the "scaling exponents" C„ and their universality with 
respect to the nature of the forcing /. This question poses 
serious difficulties since it is impossible to derive a closed 
form theory for a given order structure function Sn, since 
any such theory involves higher order unfused correlation 
functions with integrations over the time variable [3, 4]. 

A closely related question with lesser theoretical diffi- 
culties pertains to other fields that couple to the velocity 
field, with the "passive scalar" case drawing most atten- 
tion during the nineties. A passive scaler </)(r, t) is a field 
that is advected by a turbulent field which itself is unaf- 
fected by it. For example, 

^+wV</) = /cVV + /. (5) 

If u and / are homogeneous and isotropic, and Re^ oo 
and K — » 0, the structure functions Sn = {[4>{r + R) — 
(/)(r)]") are stipulated to be homogeneous fimctions of 
the their arguments with scaling exponents 

Dimensional considerations predict (n = £,n = n/S, 
with Cs = 1 being an exact result from fiuid mechanics, 
going back to Kolmogorov [5]. Experimental and simu- 
lations data deviated from these predictions (except, of 
course, for n = 3), and a hot pursuit for an example 
where these exponents could be calculated theoretically 
was inevitable. The first example that yielded to anal- 
ysis was the Kraichnan model [6], in which u is not a 
generic velocity field, but rather a random Gaussian field 
whose second order structure function scales with a scal- 
ing exponent ^2 as in Eq. (4), but is (5-correlated in time. 
This feature of the advecting field leads to a great theo- 
retical simplification, not as much as to provide a closed 
form theory for Sn, but enough to allow a derivation of 
a differential equation for the simultaneous 2n-th order 
correlation function J^2n = (?^(^i) • • • ^(i"2n)), having the 
symbolic form [6] 

OT2n = RHS{T2n-2) ■ (6) 

Guessing the scaling exponent of F2n by power counting 
and balancing the LHS against the RHS yields dimen- 
sional scaling estimates which, in this case, are ^2?i = 
(2 — C2)n. The crucial observation, however, is that the 
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differential equations (6) possess homogeneous solutions 
of the equation OT^n = [7, 8]. These "zero modes" 
are homogeneous functions of their arguments but their 
exponent cannot be guessed from power counting; the 
scaling exponents are anomalous — i.e., ^2ri < (2 — (^2)n — 
and therefore dominant at small scales. The exponents 
could be computed in perturbation theory around (^2 = 0, 
demonstrating for the first time that dimensional scaling 
exponents are not the solution to the problem. 

An appealing interpretation of the physical mechanism 
for anomalous exponents of the Kraichnan model was 
presented in the framework of the Lagrangian formula- 
tion [9]. In this formulation an nth order correlation 
function results from averaging over all the Lagrangian 
trajectories of groups of n fluid points that started some- 
where at t = — oo and ended their trajectories at points 
ri . . . r„ at time t = Q. Analyzing this dynamics it 
turned out that the Richardson diffusion of these groups 
did not contribute anything to the anomalous scaling. 
Rather, it is the dynamics of the shapes (triangles for 3 
points, tetrahedra for 4 points etc) that is responsible for 
the anomaly. In fact, the anomalous scaling exponents 
could be related to eigenvalues of operators made from 
the shape-to-shape transition probability [10]. The zero 
modes discussed above are distributions over the space 
of shapes that remain invariant to the dynamics [11]. It 
appears that these findings of the importance of shapes 
rather than scales in determining anomalous exponents 
is a new contribution to the plethora of anomalous expo- 
nents in field theory, and it would be surprising if other 
examples where shapes rather than scales are crucial will 
not pop-up in other corners of field theory, classical or 
quantum- mechanical. 

The finding of distributions that remain invariant to 

the dynamics meant that also in the Eulerian frame there 
must be such distributions, since the change from La- 
grangian to Eulerian is just a smooth change of coordi- 
nates. Indeed this was the case, and this was the clue 
how to generalize the results of the non-generic Kraich- 
nan model to the generic case Eq. (5) with a generic 
velocity field that stems from the Navier-Stokes equa- 
tions. The central comment is that the decaying passive 
scalar problem, i.e. Eq (5) with / — is a linear problem 
for which one can always define a propagator from J^n 
at t = (i.e. {(j){r\,t = 0) . . . (^(r„, t = 0)) to the same 
object at time t (note that for the decaying problem this 
is no longer a stationary quantity) [12]. This propagator 
possesses eigenfunctions of eigenvalue 1 which are homo- 
geneous functions of their arguments, characterized by 
anomalous exponents. These are the analogs of the zero 
modes of the Kraichnan model, and they are responsible 
for the anomalous exponents in the generic case [13, 14]. 
Thus the general statement that can be made is that 
the anomaly in the case of passive scalar, generic or not, 
is due to the existence of "statistically preserved struc- 
tures"; these can increases or decrease in every single 
experiment, but remain invariant on the average. This is 
a novel notion that pertains to nonlinear nonequilibrium 



systems without a known analog in equilibrium statistical 

physics. 

At present it is still unclear whether the insight gained 
from linear models might have direct relevance to the 
nonlinear problem itself. Some positive indications in 
this direction can be found in [15] , but much more needs 
to be done here before conclusions can be drawn. 



III. STATISTICAL THEORY OF ANISOTROPIC 
HOMOGENEOUS TURBULENCE 

As mentioned above, the very same agents that 
force turbulence tend to destroy also homogeneity and 
isotropy. In this section we are concerned about the loss 
of isotropy and review the extensive work that has been 
done to come to grips with this issue in a systematic fash- 
ion. Since this subject has been reviewed extensively [16], 
we limit this section to only a few essential comments. 

The need for rethinking the issue of loss of isotropy 
was underlined by the appearance of several papers where 
anisotropic flows were analyzed disregarding anisotropy, 
and exponents were extracted from data assuming that 
the inertial range scales were isotropic. The results were 
confusing: scaling exponents varied from experiment to 
experiment, and from one position in the flow to another. 
If this were indeed the case, the notion of universality in 
turbulence would fail irreversibly. In fact, it can now be 
shown that all these worrisome results can be attributed 
to anisotropic contributions in the inertial range, as ex- 
plained below. 

The basic idea in dealing with anisotropy is that the 
equations of fluid mechanics are invariant to all rotations. 
Of course, these equations are also nonlinear, and there- 
fore one cannot foliate them into the sectors of the SO (3) 
symmetry group. The equations for correlation functions 
are, however, linear (though forming an infinite hierar- 
chy). Thus by expanding the correlation functions in the 
irreducible representations of the symmetry group, one 
gets a set of equations that are valid sector by sector 
[17]. The irreducible representations of the S0(3) sym- 
metry group are organized by two quantum numbers j, m 
with j = 0, 1, 2, . . . and m = —j, —j + 1, . . .j. It turns 
out that the m components are mixed by the equations of 
motion, but the j components are not. Accordingly one 
can show that an n-point correlation function admits and 
expansion 

-Fn(n, r2, . . . r„) = ^ Agjm{ri,r2, r„)-Bgj>„(f i, fa ...?„), 

qjm 

(7) 

where f is a unit vector in the direction of r, and Agjm 
is a homogeneous function of the scalar ri . . . r„, 

^gjm(Ari,Ar2,...,Ar„) = A^"'^,j>„(ri,r2, . . . ,r„) . 

(8) 

Here Cn^ is the scaling exponent characterizing the j- 

sector of the symmetry group for the nth order correla- 
tion function. Bq^mi^i, f 2 • ■ . f„) are the n-rank tensorial 
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irreducible representations of the SO (3) symmetry group, 
and the index q in Eq. (7) is due to the fact that higher 
order tensors have more than one irreducible representa- 
tion with the same jm [17]. 

It was shown that this property of the n-th order cor- 
relation functions is inherited by the structure functions 
as well [18]. Since these are scalar functions of a vector 
argument they get expanded in standard spherical har- 



monics 



with 



.(R) 



Sn{R) ='^ajm{r)(t>jm{'R-) , 



^{r) 



(9) 



(10) 



The main issue for research was the numerical values of 
this plethora of scaling exponents. The subject was re- 
viewed in depth in [16], and we therefore limit our com- 
ments here to just the bare essentials. 

Of considerable help in organizing the scaling expo- 
nents in the various sectors of the symmetry group where 
the Kraichnan model and related models (like the passive 
vector model with pressure), where the exponents could 
be computed analytically in the Eulerian frame in any 
sector of the symmetry group. The central quantitative 
result of the Eulerian calculation is the expression for the 
scaling exponent which is associated with the scaling 
behavior of the n-order correlation function (or structure 
function) of the scalar field in the j'-th sector of the sym- 
metry group. In other words, this is the scaling exponent 
of the projection of the correlation function on the j-th 
irreducible representation of the SO(d) symmetry group, 
with n and j taking on even values only, n = 0, 2, . . . and 
J = 0,2,...' [20]: 



n{n + d) {d+l)j{j + d-2) 



2{d + 2) 2{d + 2){d-l) 



+ 0{e'). 
(11) 

The result is valid for any even j < n, and to 0(e). In 
the isotropic sector (j = 0) we recover the result of [7]. 
It is noteworthy that for higher values of j the discrete 
spectrum is a strictly increasing function of j. This is 
important, since it shows that for diminishing scales the 
higher order scaling exponents become irrelevant, and for 
sufficiently small scales only the isotropic contribution 
survives. As the scaling exponent appear in power laws 
of the type (r/A)^, with A being some typical outer scale 
and r <C A, the larger is the exponent, the faster is the de- 
cay of the contribution as the scale r diminishes. This is 
precisely how the isotropization of the small scales takes 
place, and the higher order exponents describe the rate 
of isotropization. Nevertheless for intermediate scales, or 
for finite values of the Reynolds and Peclet numbers, the 
lower lying scaling exponents will appear in all the mea- 
sured quantities, and understanding their role and disen- 
tangling the various contributions cannot be avoided. 

For Navier-Stokes turbulence the exponents cannot be 
computed analytically, but the results that are obtained 



from analyzing both experiments [18] and simulations 
[19] indicate that the picture obtained for the Kraichnan 
model repeats. The isotropic sector is always leading (in 
the sense that sealing exponents belonging to higher sec- 
tor are numerically larger). There is growing evidence 
of universality of scaling exponents in all the sectors, 
but this issue is far from being settled, and more exper- 
iments and simulations are necessary to provide decisive 
evidence. It is noteworthy that the issue of universality 
of the exponents in the isotropic sector is here expanded 
many- fold into all the sectors of the symmetry group, and 
this is certainly worth further study. 



IV. WALL-BOUNDED TURBULENCE 

Turbulent flows of highest relevance for engineering 
application possess neither isotropy nor homogeneity. 
For example, turbulent flows in channels and pipes are 
strongly anisotropic and inhomogeneous; indeed, in a sta- 
tionary plane channel flow with a constant pressure gra- 
dient p' = —dp/dx the only component of the mean ve- 
locity V, the streamwisc component Vx = V, depends 
strongly on the wall normal direction z; so do the deriva- 
tives of Vx with respect to z and the second order quan- 
tities such as mean-square-fluctuations. A long-standing 
challenge for engineers is the description of the pro- 
flies of the mean velocity and second order fluctuations 
throughout the channel or pipe at relatively high but fi- 
nite Reynolds numbers. 

To understand the issue, focus on a channel of width 
2L between its parallel walls, where the incompressible 
fluid velocity U{r,t) is decomposed into its mean (i.e., 
average over time) and a fluctuating part 

U{r, t) = V{r) + u{r, t) , V{r) = {U{r, t)) . (12) 

Near the wall, the mean velocity proflles for different 
Reynolds numbers exhibit (to the lowest order) data col- 
lapse once presented in wall units, where the Reynolds 
number Re^, the normalized distance from the wall z~^ 
and the normalized mean velocity V~^{z~^) are defined 
(for channels) by 



Re,- 



y+ — 



= zRSr/L , V+ = V/y 



The classical theory of Prandtl and von-Karman for 
infinitely large Rei- is based on dimensional reasoning 
and on the assumption that the single characteristic scale 

in the problem, is proportional to the distance from the 
(nearest) wall (and see below for details). It leads to the 
celebrated von-Karman log-law [1] 



V+iz+) = K-Hii{z+) + B, 



(13) 



which serves as a basis for the parametrization of turbu- 
lent flows near a wall in many engineering applications. 

On the face of it this law agrees with the data (see, e.g. 
Fig. 4) for relatively large z~^, say for z~^ > 100, giving 
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K ^ 0.4 and B ^ 5. The range of validly of the log-law 
is definitely restricted by the requirement ^ <C 1, where 
( = z/L (channel) oy ( = r/R (pipe of radius R). For 
^ ~ 1 the global geometry becomes important leading to 
unavoidable deviations of V^{() from the log-law (13), 
known as the wake. 

The problem is that for finite Re^ the corrections to 
the log-law (13) are in powers of e = l/lnRe,- [21] and 
definitely cannot be neglected for the currently largest 
available direct numerical simulation (DNS) of channel 
flows (Re^ = 2003 [22, 23] or £ w 0.13). Even for Re,- 
approaching 500, 000 as in the Princeton Superpipe ex- 
periment [24], e ft! 0.08. This opens a Pandora box with 
various possibilities to revise the log-law (13) and to re- 
place it, as was suggested in [21], by a power law 



V+{z+) = C{Rer){z+y(^^-^ . 



(14) 



Here both the coefficients C(Rei-) and the exponents 
7 (Re,-) were represented as asymptotic series expansions 
in e. The relative complexity of this proposition com- 
pared to the simplicity of Eq. (13) resulted in a less than 
enthusiastic response in the fluid mechanics community 
[25], leading to a rather flerce controversy between the 
log-law camp and the power-law camp. Various attempts 
[21, 24-28] to validate the log-law (13) or the alternative 
power-law (14) were based on extensive analysis of exper- 
imental data used to fit the velocity profiles as a formal 
expansion in inverse powers of s or as composite expan- 
sions in both 2;+ and (. 

Recently a complementary approach to this issue was 
proposed on the basis of experience with critical phe- 
nomena where one employs scaling functions rather than 
scaling laws [29] . The essence of this approach is the real- 
ization that a characteristic scale, say £, may depend on 
the position in the flow. The simple scaling assumption 



near the wall, i?+ 



leads to the log-law (13). The 



alternative suggestion of [21], oc (2;+)"(R'ex)^ leads to 
alternative power-law (14). But there is no physical rea- 
son why i should behave in either manner. Instead, it was 
shown that £/L should depend on ( = z/L, approaching 
k( in the limit C ~* (in accordance to the classical 
thinking). However for ( ^ 1, £ should saturate at some 
level below kL due to the effect of other walls. We re- 
call now the recent analysis of DNS data that provides a 
strong support to this idea, allowing one to get, within 
the traditional (second-order) closure procedure, a quan- 
titative description of the following three quantities: the 
mean shear, S{z) = dV{z)/dz, the kinetic energy den- 
sity (per unit mass), K{z) = (|ix|^)/2, and the tangential 
Reynolds stress, W{z) = — {uxUz}- This is achieved in 
the entire flow and in a wide region of Re,-, using only 
three Re^-independent parameters. 

The first relation between these objects follows from 
the Navier-Stokes equation for the mean velocity. The 
resulting equation is exact, being the mechanical balance 
between the momentum generated at distance z from the 
wall, i.e. p'(L — z), and the momentum transferred to the 



wall by kinematic viscosity and turbulent transport. In 
physical and wall units it has the form: 

uS + W= p'{L -z) ^ S+ + W+ = l-C, . (15) 

Neglecting the turbulent diffusion of energy (known to be 
relatively small in the log- low region), one gets a second 
relation as a local balance between the turbulent energy 
generated by the mean flow at a rate SW , and the dissi- 
pation at a rate = j^(|Vii|^): ~ SW. In stationary 
conditions e^^ equals also the energy flux from the outer 
scale of turbulence, , toward smaller scales. Thus flux 
is estimated as ^j^{z)K{z), where 7^(2;) is the typical 
eddy turn over inverse time, estimated as ^jK(z)/ij^(z). 
This gives rise to the other (now approximate) relations: 



5+^+ « £+ 



(16) 



The third required relationship can be obtained from the 

Navier Stokes equation, similarly to Eq. (16), as the local 
balance between the rate of Reynolds stress production 



SK and its dissipation 



SK. The main 



contribution to comes from the so-called Return-to- 
Isotropy process and can be estimated [30], similarly to 
£^ , as 7,^ with = \fK/l^ , involving yet another 
length-scale which is of the same order of magnitude 
as £j^. Thus one has, similarly to Eq. (16): 



S+K+^s+, e+='y+W+ = W+ 



(17) 



Now we show that the source of confusion is the assump- 
tion that the length scales can be determined a priori as 
£'^ „ (z+)" with a = 1 or a 7^ 1. In reality we have an- 
other characteristic length-scale, i.e. L, that also should 
enter the game when C = z/L is not very small. The 
actual dependence £yy and £j^ on z and L can be found 
from the data. Consider first deflned by Eq. (17), 



and introduce a new scale = £^„t^ 



(2"'")/Kjy such that 



.+ _ iy+(2+,Re.)r^(2+) 

S+{z+, Rer)^/K+{z+,Rer) 



(18) 



Here, r^{z~^) is a universal, i.e. Re^-independent di- 

mensionless function of z"*", chosen such that new scale 
£„ / L = £^ /Rot- becomes a Re^-independent function 
of only one variable ( = {z/L) = (z^/Rex)- The di- 
mensionless constant » 0.20 is chosen to ensure that 
Ymiz^L £^{C,) = z'^ . Note that if r„ were a constant, 
£„ would have started near the wall quadratically, i.e. 
as 2; X Later i'+ would have become oc 2;+ for 

50 <i; 2+ <C Re,- [30]. Thus to normalize it to slope 1 
we need the function (2;+) that behaves oc l/z"*" for 
z'^ <C 50 and approaches unity (under a proper choice 
of K^) for 0+ 50. A choice that leads to good data 
collapse reads 



r^(z+)= l+{£lJz+)' 



1/6 



t 



buf ' 



49, 



(19) 
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FIG. 1: Color online. Comparison of the theoretical mean velocity profiles (red solid lines) at different values of Re^ with the 
DNS data for the channel flow [22, 23] (Left panel, grey squares; model with l^uf = 49, k = 0.415, = 0.311) and with the 
experimental Super-Pipe data [24] (Right panel, grey circles; model with ^buf = 46, k = 0.405, is = 0.275). In orange dashed 
hne we plot the viscous solution = . In green dashed dotted line we present the von-Karman log-law. Note that the 
theoretical predictions with three Rer-independent parameters fits the data throughout the channel and pipe, from the viscous 
scale, through the buffer layer, the log-layer and the wake. For clarity the consequent plots are shifted vertically on five units. 
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FIG. 2: Color online. The scaling function ^^(C)/Rer (Left panel), £j((^)/ReT (Middle panel) and the final scaling function 
£+(C) (Right panel), as a function of = zjL, for four different values of Re^, computed from the DNS data [22, 23]. Note 
the data collapse everywhere except at ^ 1 where VK"*" ~ S"*" <C 1 and accuracy is lost. The green dash line represents 
= (1 — ^/2) with a saturation level 0.5; in orange solid line we show the fitted function Eq. (24) with ^sat = 0.311. 



where is a Rer-independent length that plays a role 
of the crossover scale (in wall units) between the buffer 
and log-law region. The quality of the data collapse for 
this scaling function is demonstrated in Fig. 2. 

The second length-scale, is determined by Eq. (16): 



(j^+(z+,Re.))3/^ 
£+(z+,ReO 



3.7 



(20) 



In Fig. 2 we demonstrate that this simple scaling function 
leads to good data collapse everywhere except maybe in 
the viscous layer. We will see below that this has only 
negligible effects on our results. 

Solution and Velocity Profiles: Solving Eqs. (16), 



(17) and accounting for Eqs. (18), (20) we find 



(21) 



where we have defined the von-Karman constant and the 
crucial scaling function {C) as follows 

^^(k3^«J1/4« 0.415, £+^[C'(C)4(C)]'/'- (22) 

The convincing data collapse for the resulting function 
^+(C)/Rer is shown in Fig. 2, rightmost panel. Substi- 
tuting Eq. (21) in Eq. (15) we find a quadratic equation 
for S with a solution: 



5+ = 



^i + (i-C)[2..^+(C)]Vv(^+)^/^ 

2[k^+(C)]V^v.(^+)'/' 



1 



(23) 
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FIG. 3: The Reynolds stress profiles (solid lines) at Rbt from 
394 to 2003 (in channel) and from 5050 to 165,000 (in pipe) in 
comparison with available DNS data (dots) for the channel. 

To integrate this equation and find the mean velocity 
profile for any value of Re,- we need to determine the 
scaling function ^+(C) from the data. A careful analysis 
of the DNS data allows us to find a good one-parameter 
fit for£+(C) [31], 

?^--{— [-1(^4)]} 

where C = C(l - C/2) and 4 ~ 0.311. The quality of the 
fit is obvious from the continuous line in the rightmost 
panel of Fig. 2. 

Finally the theory for the mean velocity contains three 
parameters, namely is together with i^^g and k. We 
demonstrate now that with these three parameters we 
can determine the mean velocity profile for any value 
RCt, throughout the channel, including the viscous layer, 
the buffer sub-layer, the log-law region and the wake. 
Examples of the integration of Eq. (23) are shown in 
Fig. 4. We trust that irrespective of the past adherence 
to the log-law camp or the power-law camp, the sympa- 
thetic reader should agree that these fits are very good. 
It remains now to estimate, using the explicit result (23), 
when do we expect to see a log-law and when the devia- 
tions due to a finite value of Rct would seem important. 
In addition, our theory results also in the kinetic energy, 
and Reynolds stress profiles which are in a qualitative 
agreement with the DNS data; for W profiles see Fig. 3. 

To show that the present approach is quite general, 
we apply it now to the experimental data that were at 
the center of the controversy [21], i.e., the Princeton Su- 
perpipe data [24]. In Fig. 4 right panel we show the 
mean velocity profiles as measured in the Superpipe com- 
pared with our prediction using the same scaling function 

(^) . Note that the data spans values of Re^ from 5050 
to 165000, and the fits with only three Rer-independent 
constants are excellent. Note the 2% difference in the 
value of K between the DNS and the experimental data; 
we do not know at this point whether this stems from 



inaccuracies in the DNS or the experimental data, or 
whether turbulent flows in different geometries have dif- 
ferent values of k. While the latter is theoretically ques- 
tionable, we cannot exclude this possibility until a better 
understanding of how to compute k from first principles 
is achieved. 

So far we discussed turbulent channel and pipe flows 
and demonstrated the existence and usefulness of a scal- 
ing function ^^(C) which allows us to get the profiles of 
the mean velocities for all values of Re,- and throughout 
the channel. While this function begins near the wall as 
z+, it saturates later, and its full functional dependence 
on C is crucial for finding the correct mean velocity pro- 
files. The approach also allows us to delineate the ac- 
curacy of the log-law presentation, which depends on z+ 
and the value of Rbt- For asymptotically large Re,- the 
region of the log-law can be very large, but nevertheless 
it breaks down near the mid channel and near the buffer 
layer, where correction to the log-law were presented. 

The future challenge is to apply this idea to other ex- 
amples of wall-bounded fiows, including time developing 
boundary layers, turbulent flows with temperature gradi- 
ents or laden with particles. There may be more typical 
"lengths" in such systems, and it is very likely that turn- 
ing these lengths into scaling functions will provide new 
insights and better models for a variety of engineering 
applications. 



V. DRAG REDUCTION BY ADDITIVES 

One severe technological problem with turbulent flows 
is that they cost a lot to maintain; the drag that the 
fluid exherts on the wall increases signiflcantly when tur- 
bulence sets in. It is therefore important that there exist 
additives, in particular polymers and bubbles, that can 
reduce this drag significantly. Over the last few years 
there had been great progress in understanding these 
phenomena, and here we provide a short review of this 
progress. 



A. Drag Reduction by Polymers 

The addition of few tens of parts per million (by 
weight) of long-chain polymers to turbulent fluid flows 
in channels or pipes can bring about a reduction of the 
friction drag by up to 80% [32-35] . This phenomenon of 
"drag reduction" is well documented and is used in tech- 
nological applications from flre engines (allowing a water 
jet to reach high floors) to oil pipes. In spite of a large 
amount of experimental and simulations data, the funda- 
mental mechanism for drag reduction has remained un- 
der debate for a long time [35-37] . In such wall-bounded 
turbulence, the drag is caused by momentum dissipation 
at the walls. For Newtonian fiows (in which the kine- 
matic viscosity is constant) the momentum flux is dom- 
inated by the so-called Reynolds stress, leading to the 
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FIG. 4: Mean normalized velocity profiles as a function of the 
normalized distance from the wall during drag reduction. The 
data points from numerical simulations (green circles) [50] and 
the experimental points (open circles) [51] represent the New- 
tonian results. The black solid line is the universal Newtonian 
line which for large agrees with von-Karman's logarithmic 
law of the wall (13). The red data points (squares) [52] repre- 
sent the Maximum Drag Reduction (MDR) asymptote. The 
dashed red curve represents our theory for the profile which 
for large y"*" agrees with the universal law (25). The blue filled 
triangles [52] and green open triangles [53] represent the cross 
over, for intermediate concentrations of the polymer, from the 
MDR asymptote to the Newtonian plug. Our theory is not 
detailed enough to capture this cross over properly. 



logarithmic (von-Karman) dependence of the mean ve- 
locity on the distance from the wall [30]. However, with 
polymers, the drag reduction entails a change in the von- 
Karman log law such that a much higher mean velocity is 
achieved. In particular, for high concentrations of poly- 
mers, a regime of maximum drag reduction is attained 
(the "MDR asymptote"), independent of the chemical 
identity of the polymer [33], see Fig. 4. During the 
last few years the fundamental mechanism for this phe- 
nomenon was elucidated: while momentum is produced 
at a fixed rate by the forcing, polymer stretching results 
in a suppression of the momentum flux from the bulk 
to the wall. Accordingly the mean velocity in the chan- 
nel must increase. It was shown that polymer stretch- 
ing results in an effective viscosity that increases linearly 
with the distance from the wall. The MDR asymptote is 
consistent with the largest possible such linear increase 
in viscosity for which turbulent solutions still exist. In 
other words, the MDR is an edge solution separating tur- 
bulent from laminar flows. This insight allowed one to 
derive the MDR as a new logarithmic law for the mean 
velocity with a slope that fits existing numerical and ex- 
perimental data. The law is universal, explaining the 
MDR asymptote. 



B. Short review of the theory 

The riddle of drag reduction can be introduced by a 
juxtaposition of the effect of polymers with respect to the 
universal law Newtonian law (13). In the presence of long 
chain polymers the mean velocity profile V^{y^) (for a 
fixed value of p' and channel geometry) changes dramat- 
ically. For sufficiently large concentration of polymers 
V^{y'^) saturates to a new (universal, polymer indepen- 
dent) "law of the wall" [33], 



V+{y+) 



In 



for y+ > 10 



(25) 



This law, which was discovered experimentally by Virk 
(and hence the notation k^), is also claimed to be univer- 
sal, independent of the Newtonian fiuid and the nature 
of the polymer additive, including flexible and rigid poly- 
mers [34]. Previous to our work in this network, the nu- 
merical value of the coefficient was known only from 
experiments, k,^~^ ~ 11.7, giving a phenomenological 
MDR law in the form [33] 



V+{y+) = 11.71ny+ - 17 . 



(26) 



For smaller concentration of polymers the situation is 
as shown in Fig. 4. The Newtonian law of the wall (13) 
is the black solid line for y"*" > 30. The MDR asymptote 
(25) is the dashed red line. For intermediate concentra- 
tions the mean velocity profile starts along the asymp- 
totic law (25), and then crosses over to the so called 
"Newtonian plug" with a Newtonian logarithmic slope 
identical to the inverse of von-Karman's constant. The 
region of values of y+ in which the asymptotic law (25) 
prevails was termed "the elastic sublayer" [33]. The rel- 
ative increase of the mean velocity (for a given p') due to 
the existence of the new law of the wall (25) is the phe- 
nomenon of drag reduction. Thus the main theoretical 
challenge is to understand the origin of the new law (25), 
and in particular its universality, or independence of the 
polymer used. A secondary challenge is to understand 
the concentration dependent cross over back to the New- 
tonian plug. In our work we argue that the phenomenon 
can be understood mainly by the influence of the polymer 
stretching on the ?/"'"-dependent effective viscosity. The 
latter becomes a crucial agent in carrying the momentum 
flux from the bulk of the channel to the walls (where the 
momentum is dissipated by friction). In the Newtonian 
case the viscosity has a negligible role in carrying the 
momentum flux; this difference gives rise to the change 
of Eq. (13) in favor of Eq. (25) which we derive below. 

The equations of motion of polymer solutions are writ- 
ten in the FENE-P approximation [38, 39] by coupling 
the fluid velocity u{r,t) to the tensor flcld of "poly- 
mer conformation tensor" R{r,t). The latter is made 
from the "end-to-end" separation vector as Rap{r,t) = 
(Taffs), and it satisfies the equation of motion 



dt 



dUg 

1 



[P(r, i)i?„0 - pIS^p] 
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P{r,t) = {pl,-pl)/{pl,-R^^ 



(27) 



and Po refer to the maximal and the equilibrium val- 
ues of the trace R^-y. In most applications pm ^ Po 

P{r,t)^ (1/(1 - ai^^T,) 

where a = The equation for the fluid velocity field 
gains a new stress tensor: 

^ + {UjV^)Ua = -VaP + Z^sV^Ua + V-yT^-y (28) 



The cocfScicnts and are bounded from above by 
unity. (The proof is |c| = \W\/K < 2\{u^Uy)\/ {ul+uD < 
1, because {u^ ± Wy)^ > 0.) 

To proceed, one needs to estimate the various com- 
ponents of the polymer conformation tensor. This was 
done in [41] with the final result that for Cp large (where 
P « 1), and Deborah number De = TS{y) ^ 1 the con- 
formation tensor is highly anisotropic, 



R{y) ~ Ryy{y) 



2Be\y) 
Be{y) 




De(y) 
1 







C{y) 



P{r,t) 
Pi 



Ra/3{r,t) - Sal3 



(29) 



Here Vg is the viscosity of the neat fluid, and Vp is a 
viscosity parameter which is related to the concentration 
of the polymer, i.e. Vpjvs ~ Cp. 

We shall use the approximation 



V, 



ta/3 



o -ft, 

^ pI 



■a/3- 



Armed with the equation for the viscoelastic medium 
we establish the mechanism of drag reduction following 
the standard strategy of Reynolds. Eq. (15) changes now 
to another exact relation [40] between the objects S and 
W which includes the effect of the polymers: 



W^vS^ '^{PR^y){y) = p'{L - y) 



(30) 



On the RHS of this equation we see the production of mo- 
mentum flux due to the pressure gradient; on the LHS we 
have the Reynolds stress, the Newtonian viscous contri- 
bution to the momentum fltix, and the polymer contribu- 
tion to the momentum flux. A second relation between 
S{y), W{y), K{y) and R{y) is obtained from the energy 
balance. In Newtonian fluids the energy is created by 
the large scale motions at a rate of W{y)S{y). It is cas- 
caded down the scales by a flux of energy, and is finally 
dissipated at a rate e, where e = j/(|Vm|^). In viscoelas- 
tic flows one has an additional contribution due to the 
polymers. Our calculation [40] showed that the energy 
balance equation takes the form: 



K , K^l^ 
av-^r+b 

y 



+ ^{Pf{{Ryy) + {R..)) = WS. (31) 

We note that contrary to Eq. (30) which is exact, Eq.(31) 
is not exact. We expect it however to yield good order 

of magnitude estimates as is demonstrated below. Fi- 
nally, we quote the experimental fact [33] that outside 
the viscous boundary layer 



W{y) 
K{y) 



c? , for Newtonian flow, 



, for viscoelastic flow. 



(32) 



The important conclusion is that the term proportional 
to {Ryy) in Eq. (31) can be written as Vp{TZyy){y)S{y). 
Defining the "effective viscosity" jy(y) according to 

i^{y) =1^0 + ^p{T^vv){y) , (33) 
The momentum balance equation attains the form 



iyiy)Siy) + W{y)=p'L . 



(34) 



It was shown in [40] that also the energy balance equation 
can be rewritten with the very same effective viscosity, 



I.e., 



i^iy) {^yK{y) + b/EMxiy) = W{y)S{y) . (35) 

In the MDR region the first term on the RHS of in Eqs. 
(34) and (35) dominate; from the first equation u{y) ~ 
1/S{y). Put in Eq. (35) this leads to S{y) ^ 1/y, which 
translates to the new logarithmic law which is the MDR. 
We will determine the actual slope momentarily. At this 
point one needs to stress that this results means that 
nu{y) must be proportional to y in the MDR regime. 
This linear dependence of the effective viscosity is one 
of the central discoveries of our approach. Translated 
back, it predicts that {Ryy) ~ y outside the boundary 
layer. This prediction is well supported by numerical 
simulations. 

The crucial new insight that explained the universal- 
ity of the MDR and furnished the basis for its calcula- 
tion is that the MDR is a marginal flow state of wall- 
bounded turbulence: attempting to increase S{y)) be- 
yond the MDR results in the collapse of the turbulent 
solutions in favor of a stable laminar solution W — 0. As 
such, the MDR is universal by definition, and the only 
question is whether a polymeric (or other additive) can 
supply the particular effective viscosity i'{y) that drives 
Eqs. (34) and (35) to attain the marginal solution that 
maximizes the velocity profile. We predict that the same 
marginal state will exist in numerical solutions of the 
Navier-Stokes equations furnished with a y-dcpendcnt 
viscosity i^iy). There will be no turbulent solutions with 
velocity profiles higher than the MDR. 

To sec this explicitly, we first rewrite the balance equa- 
tions in wall units. For constant viscosity (i.e. v{y) = 
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vq), Eqs. (34)-(35) form a closed set of equations for 
S+ = Svo/{P'L) and W+ = W/VWZ in terms of two 
dimensionless constant i5+ = a^jK/W (the thickness of 
the viscous boundary layer) and = 6/cy (the von 
Karman constant). Newtonian experiments and simula- 
tions agree well with a fit using 5+ ~ 6 and ^ 0.436 
(see the black continuos line in Fig. 4 which shows the 
mean velocity profile using these very constants). Once 
the effective viscosity vd/) is no longer constant we ex- 
pect cv to change and consequently the two dimension- 
less constants will change as well. We will denote the new 
constants as A and respectively. Clearly one must re- 
quire that for y{y)/yQ ^ 1, A ^ 5+ and — » k^. The 
balance equations are now written as 



^ iv )— + 



W+{y+) = 1 
S+ . 



(36) 
(37) 



where v^{y'^) = v{y^)/vQ. Substituting now 5+ from 
Eq. (36) into Eq. (37) leads to a quadratic equation for 
\/W^ . This equation has as a zero solution for W+ (lam- 
inar solution) as long as i^+(y^)A/y"'" = 1. Turbulent 
solutions are possible only when ;/+(?/+)A/y"'" < 1. Thus 
at the edge of existence of turbulent solutions we find 
oc y^ for y+ 1. This is not surprising, since it was 
observed already in previous work that the MDR solution 
is consistent with an effective viscosity which is asymp- 
totically linear in y"*" [42, 43]. It is therefore sufficient to 
seek the edge solution of the velocity profile with respect 
to linear viscosity profiles, and we rewrite Eqs. (36) and 
(37) with an effective viscosity that depends linearly on 
y+ outside the boundary layer of thickness 



[1 + a{y+ - S+)]S+ + W+ = 1, 



(38) 

(39) 



We now endow A with an explicit dependence on the 
slope of the effective viscosity y^{y), A = A (a). Since 
drag reduction must involve a decrease in W^, we expect 
the ratio a'^K/W to depend on a, with the constraint 
that A(a) S+ when a ^ 0. Although A, 6+ and 
a are all dimensionless quantities, physically A and 5+ 
represent (viscous) length scales (for the linear viscos- 
ity profile and for the Newtonian case respectively) while 
is the scale associated to the slope of the linear vis- 
cosity profile. It follows that aS~^ is dimensionless even in 
the original physical units. It is thus natural to present 
A(a) in terms of a dimensionless scaling function f{x), 



A(a) = 5+f{a5+) . 



(40) 



Obviously, /(O) = 1. In [44] it was shown that the bal- 
ance equation (38) and (39) (with the prescribed form of 

the effective viscosity profile) have an non-trivial symme- 
try that leaves them invariant under rescaling of the wall 
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FIG. 5: The solution for 10^TF+ (dashed line) and J/+5+ 
(solid line) in the asymptotic region j/^ (5+, as a function 
of a. The vertical solid line a = 1/25"*" = 1/12 which is the 
edge of turbulent solutions; Since VW+ changes sign here, 
to the right of this hne there are only laminar states. The 
horizontal solid line indicates the highest attainable value of 
the slope of the MDR logarithmic law l/zty = 12. 



units. This symmetry dictates the function A (a) in the 
form 



A(a) = 



S+ 



1-0!^+ 



(41) 



Armed with this knowledge we can now find the maximal 
possible velocity far away from the wall, y+ ;s> 5+. There 
the balance equations simplify to 



ay+S+ + T4^+ = 1 , 
aA^ia) + Vw+/kc = y+S+ 



(42) 
(43) 



T hese equations have the ?/+-independent solution for 
VW+ and y+S+: 



-S+ = aA'^{a) + Vw+/K^ 



(44) 



By using equation (44) (see Fig. 5), we obtain that the 
edge solution {W~^ 0) corresponds to the supremum 
of y~^ , which happens precisely when a = 1/A(a). 
Using Eq. (41) we find the solution a = Um = 1/2(5+. 
Then y~^S~^ = A{am), giving = 2^+. Using the 
estimate (5+ « 6 we get the final prediction for the MDR. 
Using Eq. (25) with = 12, we get 



V+{y+) « 121ny+ - 17.8 



(45) 



This result is in close agreement with the empirical law 
(26) proposed by Virk. The value of the intercept on the 
RHS of Eq. (45) follows from Eq. (25) which is based 
on matching the viscous solution to the MDR log-law 
in [42]. Note that the numbers appearing in Virk's law 
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correspond to S'^ = 5.85, which is well within the error 
bar on the value of this Newtonian parameter. Note that 
we can easily predict where the asymptotic law turns into 
the viscous layer upon the approach to the wall. We can 
consider an infinitesimal and solve Eqs. (36) and (37) 
for 5+ and the viscosity profile. The result, as before, is 
'^^{y) = ^{ctm)y~^- Since the effective viscosity cannot 
fall bellow the Newtonian limit zv+ = 1 we see that the 
MDR cannot go below y+ = A{am) = 25+. We thus 
expect an extension of the viscous layer by a factor of 2, 
in very good agreement with the experimental data. 

C. Non-universal aspects of drag reduction by 
polymers 

When the concentration of polymers is not large 
enough, or when the Reynolds number is too low, the 
MDR is attained only up to some value of y+ that de- 
pends in a non-universal manner on the Reynolds num- 
ber and on the nature of the polymer [45]. These non- 
universal turn-backs to the so called "Newtonian plug" 
can be understood theoretically, and we refer the reader 
to [46, 47] for further details. 

D. Drag reduction by micro-bubbles 

Finally, we should mention that drag reduction by 
polymers is not the solution for many technologically 
pressing problems, the most prominent of which is the 
locomotion of ships. Here a more promising possibil- 
ity is the drag reduction by bubbles, a subject that is 
much less developed than drag reduction by polymers. 
For some recent papers on this subject see for example 
[48] and references therein; we we stress that this sub- 
ject is far from being exhausted by these papers, and we 
expect that more work should appear on this subject in 
the near future. 



VI. THERMAL CONVECTION 

Convection often occurs in conjunction with other fea- 
tures such as rotation, magnetic field and particulate 
matter, so the knowledge of the subject is relevant to sev- 
eral closely related fields. The complexity of the under- 
lying equations has precluded much analytical progress, 
and the demands of computing power are such that rou- 
tine simulations of large turbulent flows has not yet been 
possible. Thus, the progress in the field has depended 
more on input from experiment, which itself has severe 
limitations. The progress in the subject, such as it is, 
has been possible only through strong interactions among 
theory, experiment and simulation. This is as it should 
be. 

The paradigm for thermal convection is the Rayleigh- 
Benard problem in which a thin fluid layer of infinite 



lateral extent is contained between two isothermal sur- 
faces with the bottom surface maintained slightly hotter. 
When the expansion coefficient is positive (as is the case 
usually), an instability develops because the hot fluid 
from below rises to the top and the colder fluid from 
above sinks to the bottom. The applied driving force is 
measured in terms of a Rayleigh number, Ra, 

Ra = gaATH^/uK, (46) 

which emerges in front of the buoyancy term and is a non- 
dimensional measure of the imposed temperature differ- 
ence across the fluid layer. Here, g is the acceleration 
due to gravity, H is the vertical distance between the 
top and bottom plates, a, v and k are, respectively, the 
isobaric thermal expansion coefficient, the kinematic vis- 
cosity and the thermal diffusivity of the fluid. Physically, 
the Rayleigh number measures the ratio of the rate of po- 
tential energy release due to buoyancy to the rate of its 
dissipation due to thermal and viscous diffusion. 

The second important parameter is the Prandtl num- 
ber 

Pr = v/k, (47) 

which is the ratio of time scales due to thermal diffusion 
{tq = I k) and momentum diffusion (r^, = H"^ /v), and 
determines the ratio of viscous and thermal boundary 
layers on the solid surfaces. With increasing Ra the dy- 
namical state of the Rayleigh-Bcnard system goes from a 
uniform and parallel roll pattern at the onset [Ra ~ 10^) 
to turbulent state at Ra ^ lO'' — 10^. (The onset value 
is independent of Pr but the latter depends strongly on 
it.) 

For purposes of theoretical simplification, it is custom- 
ary to assume that the thermal driving does not affect the 
pressure or the incompressibility condition, and that its 
only effect is to introduce buoyancy. This is the Boussi- 
nesq approximation. How closely the theoretical results 
correspond to observations depends on how closely the 
experiments obey the Boussinesq approximation. It is 
also not clear if small deviations from the ideal bound- 
ary conditions produce only small effects. 

A. Experiments using cryogenic helium 

Since many examples of convection occur at very high 
Rayleigh numbers [54] , it is of interest to understand the 
heat transport characteristics in that limit. It is also 
necessary to be able to cover a large range of Ra in or- 
der to be able to discover the applicable scaling laws. 
Cryogenic helium has been used successfully for the pur- 
pose. Though experiments in conventional fluids have 
been valuable [55, 56], the Rayleigh number has been 
pushed to the limit only through cryogenic helium. The 
same properties that make it a suitable fluid for convec- 
tion studies also makes it suitable for creating flows with 
very high Reynolds numbers [57]. 
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Historically, a small "supcrfluid wind tunnel" was con- 
structed [58] with the idea of exploiting the supcrfluid 
properties of helium II for obtaining very high Reynolds 
numbers. Potential flow was observed for low velocities, 
with no measurable lift on a pair of fly wings hanging 
in the tunnel, but the inevitable appearance of quantized 
vortices (see section VII on superfluid turbulence) altered 
that picture for higher flow speeds. Threlfall [59] recog- 
nized the advantages of using low temperature helium gas 
to study high- i?a convection. The later work by Libch- 
aber and co-workers [60] brought a broader awareness of 
the potential of helium. The work of Refs. [61, 62] is a 
natural culmination of this cumulative effort. 



Fluid 


T{K) 


P{Bar) 


a/vK 


Air 


293 


1 


0.12 


Water 


293 


1 


14 


Helium I 


2.2 


SVP 


2.3 X 10^ 


Helium II 


1.8 


SVP 




Helium gas 


5.25 


2.36 


6 X 10^ 


Helium gas 


4.4 


2 X 10"* 


6 X 10-3 



Table 1. Values of the combination of fluid properties 
a/ UK for air, water and helium. From Ref. [63]. 

The specific advantage of using helium for convection 
is the huge value of the combination a/vn near the criti- 
cal point. This can generate large Ra (see Table 1). For a 
fluid layer some 10 meters tall and a reasonable temper- 
ature difference of Q.hK, Rayleigh numbers of the order 
10^^ are possible. Table 1 also shows that a/un is quite 
small at pressures and temperatures sufficiently far away 
from the critical value. In fact, the range shown in the 
table covers a factor of 10^^, so any experiment of fixed 
size H can yield at least 12 decades of the control param- 
eter Ra by this means alone. However, if H is chosen to 
be large enough, this entire range of Ra can be shifted to 
a regime of developed turbulence where well-articulated 
scaling relations may be observed. This tunability is es- 
sentially impossible for air and water, especially because 
one cannot use more than modest temperature difference 
to increase Ra (due to the attendant non-Boussinesq ef- 
fects. Sect. VI C). For other advantages in using helium, 
see [63]. 



B. The scaling of the heat transport 

The heat transport in convection is usually given in 
terms of the Nusselt number Nu 



NU : 



Qcond 



(48) 



where q is the total heat flux, qcond is the heat flux in the 
absence of convection, given by Fourier's law, and fc/ is 
the thermal conductivity of the fluid. Nu represents the 



ratio of the effective turbulent thermal conductivity of 
the fluid to its molecular value. One goal of convection 
research is to determine the functional relation Nu = 
f{Ra,Pr). This relation is at least as fundamental as 
the skin friction relation in isothermal flows. 

Figure 6, reproduced from Ref. [63], illustrates the 
enormous range of Ra and Nu that is possible in low 
temperature experiments of modest physical size. The 
Nusselt numbers have been corrected here for sidewall 
conduction and also for finite thermal conductivity of the 
plates (and both corrections are small, see [63]). That 
one can reach Nusselt numbers as high as 10^ bears tes- 
timony to the great importance of turbulence as a subject 
of serious study. 

We have shown this figure in part because it represents 
the highest Ra achieved so far under laboratory condi- 
tions and also the largest range of Ra in the turbulent 
scaling regime, both of which represent the fulfilment of 
the promise of cryogenic helium gas. The average slope 
over 11 decades is 0.32, close to 1/3. In part, we show 
the graph because one might have hoped that such an 
unusual figure spanning many decades in Ra might have 
a finality to it. Perhaps it does. However, experiments of 
Chavanne et al. [62], and by Niemela & Sreenivasan [65] 
for a different aspect ratio, have found a scaling exponent 
rising beyond 1/3 towards the very highest Ra. The plau- 
sible conclusions of Niemela & Sreenivasan [65, 66] were 
that those data corresponded to large departures from 
Boussinesq conditions and to variable Prandtl number 
(remembering that increasing Prandtl number serves to 
stabilize and laminarize the boundary layers, in contrast 
to conditions prescribed for observing the 1/2-power scal- 
ing), but it is important to test these plausible conclu- 
sions directly. We shall momentarily discuss the current 
work in this direction. If we ignore the apparently non- 
Boussinesq regime, it has been argued in Refs. [65, 66] 
that the scaling exponent from existing data is most likely 
consistent with a value close to 1/3. 

As already mentioned, computations have not yet ap- 
proached experiments in terms of high Ra, but their ad- 
vantage is that Pr can be held constant and the Boussi- 
nesq approximation can be enforced strictly. The limit of 
computational ability has recently been pushed by Am- 
ati et al. [64], who have reached Rayleigh numbers of 
2 X 10^^. Even though this number is still about three 
orders of magnitude lower than the highest experimen- 
tal value, it has become quite competitive with respect 
to many other experiments. This work suggests that the 
one-third exponent is quite likely, reinforcing the conclu- 
sion of Refs. [65, 66]. Computational simulations have 
also explored the effects of finite conductivity, sidewall 
conduction and non-Boussinesq effects [67, 68]. 

In spite of the limitations of Ra attainable in simula- 
tions, much of the detail we know about boundary layers 
and fiuctuations come from them. Direct knowledge of 
the velocity is most desirable in understanding the dy- 
namics of plumes and boundary layers, and also the im- 
portance of the mean wind. Experiments in convection 
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have limited themselves to measuring the mean wind and 
temperature at a few points, but not the spatial struc- 
tures. The conventional techniques of velocity measure- 
ments and flow visualization are fraught with difficulties, 
as has been discussed in [63] . 

We should now discuss the contributions of the theory 
to the heat transport problem. Two limiting cases for the 
scaling of Nu have been considered. The first scenario 
imagines that the global flux of heat is determined by 
processes occurring in the two thermal boundary layers 
at the top and bottom of the heated fluid layer. Then 
the intervening turbulent fluid, being fully turbulent and 
"randomized" , acts as a thermal short circuit and there- 
fore its precise nature is immaterial to the heat flux. We 
can then determine the relation to be Nu ^ Ra^^^ [69]. 
This scaling assumes that the heat flux has no depen- 
dence on H. In the limit in which molecular properties 
are deemed irrelevant in determining heat transport — 
that is, when boundary layers cease to exist — an expo- 
nent of 1/2 (modiflo logarithmic corrections) has been 
worked out phenomenologically [70]. There has been an 
alternative theory [60] that obtains the 2/7-ths scaling 
through intermediate asymptotics, but the experimental 
result that motivated the work has not been sustained by 
more recent work. 

The upperbound theory, though quite old (see Refs. 
[69, 71]), has been taken to new levels through the efforts 
of Constantin and Doering (e.g., Ref. [72]), as well as by 
others more recently, and has contributed some valuable 
hints on the heat transport law. The latest summary is 
as follows: 

1. Arbitrary Prandtl number: Nu < Ra}/^ uniformly 
in Prandtl number [72] . This result rules out the Prandtl 

number dependence such as Pr^/"^ [73, 74] and Pr"^/"^ 
[70]. In particular, the latter paper was written when the 
boundary layer structure was understood much less, and 
there is a need to reconstruct its arguments afresh, in 
particular for the reassessment of the Rayleigh number 
at which the so-called "ultimate regime" is supposed to 
prevail for Prandtl numbers of order unity. 

2. Large but finite Prandtl number: The largeness 
of the Prandtl number is prescribed by the condition 
Pr > cRa, where c is a constant of the order unity. Un- 
der this condition, the upperbound is given by Nu < 
Ra^/^{lnRa)'^/^ [75]. For higher Rayleigh numbers the 
upperbound is still given by (1) above. 

3. Infinite Prandtl number: The latest result due to 
Doering et al. [76], is Nu < CRa^/^{lnRaf/^ . Robust 
calculations by Ireley et al. [77], which still seem to fall 
short of proof, is Nu < aRa}/^, where a is a constant of 
the order unity. 

Thus, as far as the upperbound theory goes, the Ra^/"^ 
result is permissible for Prandtl numbers of the order 
unity, though some semi-analytical results on Prandtl 
number dependencies are ruled out as noted above. 

Finally, we mention the effect of rough surfaces on the 
global heat transfer rate [78, 79] and the presence of a 
weakly organized mean wind [80-83] . These studies have 
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FIG. 6: Log-log plot of the Nusselt number versus Rayleigh 
number. The line through the data is a least-square fit over 
the entire Ra range, and represents a dlogNu/dlogRa slope 
of 0.32. 



added to our understanding of turbulent convection. The 
wind phenomenon has had a rather broad reach; e.g., 
quantitative observations of occasional reversals of the 
mean wind flow direction have been shown to be related 
to simple models of self-organized criticality [84]. Fur- 
thermore, the lifetimes of the metastable states of the 
bi-directional mean flow have intriguing analogies with 
reversals of the Earth's magnetic field polarity, a phenom- 
ena arising from turbulent convection within the outer 
core [85]; there is also a quantitative statistical analogy 
with the lifetime of solar flare activity driven by turbu- 
lent convection in the Sun's outer layer [86]. This latter 
conclusion may indicate the existence of an underlying 
universality class, or a more direct physical similarity in 
the convective processes that lead to reshuffling of the 
magnetic footprints and ultimately to flare extinction. 

C. Non-Boussinesq effects 

One possible measure of Boussinesq conditions is that 
the fractional change in density across the layer, 

^ = aAT, (49) 
P 

must be small. On the basis of a comparison to the 
Boussinesq problem at the onset of convection, it is gen- 
erally assumed that values of aAT < 0.2, or a 20% vari- 
ation of density across the flow thickness, is acceptable. 
This criterion is indeed satisfied up to very high values of 
Ra (above 10'^''' for one set of data [65] and above 10^^ for 
another [61]), although there is no assurance that asym- 
metries of this magnitude are irrelevant at such high Ra. 
In fact, a more stringent requirement by a factor of 4 was 
adopted in Ref. [65]. 
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Because of the importance of the non-Boussincsq ef- 
fects, as discussed in Ref. [65], recent attention has been 
focused on them. The earhest exploration was by Wu 
& Libchaber [87], who reported top-bottom asymmetry 
in boundary layers as a main characteristic and a drop 
in the ratio of temperature drop across top to bottom 
boundary layer as the Rayleigh number increases. Veloc- 
ity profiles measured in a follow-up paper [88], at lower 
Ra, using glycerol, also showed an asymmetry. Ahlers 
and collaborators [89] showed that the non-Boussinesq 
effects depend on the fluid, as one could expect. For wa- 
ter, Nu showed a modest decrease with increase in AT. 
For ethane, they found larger Nu than in the Boussinesq 
case, nearly 10% higher when aAT = 0.2. 

Because there are many possible non-Boussinesq ef- 
fects and their relative importance depends on the fluid 
and the oper ating conditions, it is difficult to study these 
effects systematically in experiments. A numerical com- 
putation by [90] in two dimensions, with glycerol as work- 
ing fluid, showed that effects on Nu were marginal, with 
some decrease in Nu with aAT for Ra > 10^. In [68], 
these effects have been explored in three dimensional con- 
vection, also computationally. The finding is that at 
least for conditions corresponding to cryogenic helium 
gas at modest Rayleigh numbers — while viscosity plays 
an important role in diminishing the movement of plumes 
to the interior of convection it is the coefficient of thermal 
expansion that affects heat transport most. 



D. Whither helium experiments? 

While thermal convection has been studied for quite 
some time, the recent surge of interest has been triggered 
by helium experiments- -even in theory and simulations. 
Indeed, experiments were ahead of theory and simula- 
tions about two decades ago. Since then, theory has been 
making its presence felt slowly and simulations have been 
making considerable inroads. Experiments have surely 
extended the parameter ranges, but, just as surely, they 
have not kept up the pace of sophistication. A major step 
in the understanding of the problem will occur only if a 
major improvement in experimental sophistication takes 
place. It is therefore useful to take stock of the situation 
briefly. It is perhaps useful even to raise the question 
as to whether the promise of helium is realizable in its 
entirety anytime soon. 

It has been recognized abundantly that the problem 
is with instrumentation and with probes of the desired 
temporal and spatial resolution. It is not clear to us 
that smaller probes based on the principles of standard 
thermal anemometry are the solution to the problem, 
part of which arises because the use of helium raises the 
Reynolds number of the probe itself to a higher value 
than in conventional fluids, leading to unfavorable (and 
poorly analyzed) heat transfer characteristics. 

In thermal convection flows, where some direct knowl- 
edge of the velocity would be most desirable even at scales 



much larger than the Kolmogorov length, the use of hot 
and cold wires is further complicated by the fact that 
they require a steady flow — and the mean wind is effec- 
tive only near the boundaries and also becomes weaker 
with Ra. Complications arise from the simultaneous 
presence of temperature fluctuations in the cell and the 
temperature fluctuations of the wire due the velocity fluc- 
tuation that is intended to be measured. 

Even if single-point measurements were possible suc- 
cessfully, the need to measure the entire velocity field in a 
turbulent flow remains to be addressed. While a number 
of hot wires at several points can be used to obtain some 
spatial information, there is a limit to this procedure. 
In principle, we may use the Particle Image Velocimetry 
(PIV) to obtain an entire two-dimensional section of the 
turbulent flow field at a given instant in time. In fact, 
PIV has been applied recently to liquid helium grid tur- 
bulence at 4.2K [91, 92], in counterflow turbulence [93] 
and in helium II turbulence [94]. It should be pointed 
out that because of the enormous amount of information 
in any one image, it is difficult to process very high data 
rates from this type of measurement. Consequently, time 
evolution of the flow cannot be easily obtained. 

Particle selection and injection remain the fundamen- 
tal hurdle for PIV measurements at low temperatures. 
Liquid helium has a relatively low density, and this makes 
it harder to flnd suitably buoyant particles that are also 
not too large. The use of hydrogen particles that match 
the density of helium has been the most promising step 
in this direction [94], but better control of the particle 
generation is needed to render the technique routinely 
usable. It is equally important to better understand the 
interaction of particles with the mixture of normal and 
super fluids [95]. 

The seeding of helium gas for thermal convection ex- 
periments is probably even more difficult owing to the 
large variation of the density, and its nominally small 
value, which at best is less than half that of the liquid 
phase. However, the liquid, as mentioned above, can also 
be used to attain high Ra, though at the expense of a 
large range. The compensation is that we know that the 
liquid flow can be seeded to some level of adequacy. 

Flow visualization can focus experimental and even 
theoretical — efforts, and yet this domain has not been 
well developed for cryogenic helium. We believe that 
there is a huge pay-off because most existing flow visual- 
izations in water and other room-temperature fluids are 
at low to moderate Rayleigh numbers, and the intuition 
that one derives from low Ra cannot easily be extended 
to high Ra. There are no technological barriers to per- 
fecting the present efforts — only one of integrating var- 
ious components together. We may also remark that it 
is not easy to test new particles in the actual low tem- 
perature environment. In experimental phase. White et 
al. [91] had resorted to testing in a pressurized SFq en- 
vironment, where the density could be matched to that 
of liquid helium. 

Where density gradients exist in the flow, visualiza- 
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tion can occur in the absence of tracer particles, using 
shadowgraphs (which depends on the density gradient) or 
schUeren technique (which depends on the second deriva- 
tive of the density). It has been demonstrated [96] that 
shadowgraphy can be used in hehum I to visuaHze even 
weak flows near the convective onset. A hght beam re- 
flected from the cell displays intensity variations resulting 
from convergence or divergence as a result of gradients 
in the refractive index. In the case of thermal convec- 
tion, these indicate the average temperature field. Note 
that the technique does not give local information, but 
can be used to visualize only global flows. In the case of 
large apparatus, installing an optically transparent but 
thermally conducting plates is a non-trivial task. 

For the case of turbulence under isothermal conditions, 
it would be possible to use helium 3 as a "dye marker" 
for shadowgraphs. 

Scattering of ultrasound is another method that can in 
principle be used for velocity measurements in helium. It 
can be used in the gas phase which makes it a plausible 
candidate for cryogenic convection experiments. How- 
ever, there would be substantial problems with achiev- 
ing sufficiently high signal-to-noise ratio resulting from 
a mismatch of acoustic impedance between the sound 
transducers and the helium. The work in this direction 
[97; 98] has not yet been adopted in cryogenic helium. 

In summary, one part of the promise of helium (namely 
large values and ranges of the control parameters) has 
been amply established; flows with huge values of Ra and 
Re have indeed been generated in laboratory-sized appa- 
ratus. However, the second part of the promise (of being 
able to develop versatile techniques for precise measure- 
ments of velocity and vorticity) has lagged behind sub- 
stantially, despite some impressive efforts. This is the 
aspect that needs substantial investment. 

A speculative possibility in this direction, cited in [63], 
is the laser-induced fluorescence of metastable helium 
molecules, which can be the "tracer" particles in PIV 
measurements. This could be a sensitive tool with good 
spatial resolution (only molecules in the intersecting re- 
gion of two crossed lasers would be excited to metastable 
states), and it is also possible that a single molecule could 
be detected spcctroscopically. Above IK, the molecules 
would likely move with the normal fluid component, al- 
lowing its velocity profile to be detected, provided that 
the molecules are not trapped on vortex lines. Such a 
technique would be complementary to more typical sec- 
ond sound attenuation measurements. 

Once the instrumentation issues are clearer, we need to 
seriously consider an experiment that can combine mod- 
erate aspect ratio (say, 4) with high Ra, constant Pr, and 
Boussinesq conditions. Such an experiment is probably 
not without considerable technical difficulties. A large 
scale low temperature apparatus could be constructed, 
say at a facility like CERN or BNL, where there is an 
adequate refrigeration capacity. Having a horizontal di- 
mension of, say, 5 meters or more would probably require 
some type of segmentation of the plates with multiplexing 



of the heating and temperature control. Fundamentally, 
this is no more complicated than the mirror arrays used 
in astrophysical observation. The bottom plate, which 
has a constant heat flux condition imposed, can be arbi- 
trarily thick since it can be supported from below. The 
top, temperature controlled plate would probably have 
some limitations in this respect. Estimates have been 
made of the cooling power required for a cell that is 5 m 
in diameter but also 10 m tall, and it is around 200W or 
less, which is not a severe requirement. 



VII. SUPERFLUID TURBULENCE 

We now review aspects of liquid helium below the 
lambda point, called helium II. At low velocities, he- 
lium II flows without friction but the situation changes 
when velocity exceeds a critical value: the fluid enters 
a state in which quantized vortices are formed sponta- 
neously into a self-sustained tangle — except under con- 
trolled conditions such as rotation, for which the vortices 
are all aligned with the direction of rotation. The vor- 
tex lines move about in the background of freely moving 
elementary excitations, which one may regard as some- 
thing of background "gas particles" . The vortices scat- 
ter the excitations as long as there is relative velocity 
between them and the "particles", thus generating fric- 
tion. The vortices are estimated to be of the order of an 
angstrom in diameter, and, as was recognized by Onsager 
in 1949, quantum mechanics constrains the circulation 
around them to be riK/m, where k is the Planck's con- 
stant and m is the mass of the helium atom; the integer 
n = 1 normally. However, the irrotational flow further 
away from the core of the vortices is classical. The mo- 
tion produced by a vortex tangle can be quite complex 
because of its complex geometry and is called superfluid 
turbulence. 



A. Analogy to classical turbulence 

One of the recent findings [99] is that the superfluid 
turbulence has the Kolmogorov form for the spectral 
density with a well-defined —5/3 power, independent of 
whether the fraction of the normal fluid (corresponding 
to excitations) is negligible or dominant. This result may 
not seem surprising if one takes the stand that any non- 
linearly interacting system of many scales will behave 
similarly to the classical Kolmogorov turbulence in the 
inertial range [1]. What is needed are the mechanisms 
of excitation at some large scale and dissipation at the 
small scale, with no further detail mattering in the iner- 
tial range. However, several problems come to the fore 
upon closer scrutiny. 

First the dissipation mechanism: It is generally ac- 
cepted that the short wavelength Kelvin waves are re- 
sponsible for dissipation [100]. These waves are created 
presumably by the impulse associated with the recon- 
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nection of vortices. For temperatures of 1 K and above, 
Kelvin waves arc damped out by the background exci- 
tations (the normal fluid) thus providing the dissipation 
mechanism. For lower temperatures, for which the nor- 
mal fluid is negligible, the energy is radiated away as 
sound at sufficiently small wavelengths. For radiation to 
be effective, one needs high velocities and short wave- 
lengths: modest motion of vortices will not do. Higher 
velocities are possible very close to the vortex core be- 
cause of the inverse power law of the potential veloc- 
ity field — and also because of reconnection events, which 
produce cusp-like local structures with sharply repelling 
velocities. 

Regarding the forcing scale, in experiments with a pull 
through grid in helium II [101], it is conceivable that the 
forcing is produced very similarly to that in classical tur- 
bulence, and is related to the mesh length and the time 
of evolution of the turbulence. In some simulations, the 
forcing scale cannot be defined unambiguously. For in- 
stance, in the important foray into superfluid turbulence 
that was made by Schwarz [102], it appears that forc- 
ing scale was the size of the computational box, as in 
the case of the simulations of the Taylor-Green problem 
by Nore et al. [103] and Araki et al. [104]. However, it 
also appears that reconnections play an important role 
in determining this scale. 

As another perspective on the same issue, the occur- 
rence of the —5/3 spectrum in superfluid turbulence may 
be regarded as surprising if one takes the stand that the 
key mechanism for energy transfer across scales in hydro- 
dynamic turbulence, namely vortex stretching, is absent 
in superfluid turbulence: no intensification and break- 
up due to vortex stretching is possible. It is the vor- 
tex break-up due to reconnections, not vortex stretching, 
that appears to be the key to the spectral distribution 
here. If this is true, it is an interesting to speculate about 
the central importance attached to vortex stretching in 
classical turbulence. 

To be sure, one should look closely at the veracity of 
claims about the —5/3 power law. Our view is that the 
available evidence is too fragile to sustain the claim on 
the existence of the —5/3 spectrum in experiment or sim- 
ulations. In experiments, the only real piece of evidence 
comes from Ref. [99], but at least to us it is not exactly 
clear what is being measured at the low end of the tem- 
perature (below 1 K), despite a nice assessment in Ref. 
[105]. In simulations of superfluid turbulence, the result 
is unconvincing because the computational box size is 
still a long way from reaching the size characteristic of 
those in classical turbulence. We make a strong case for 
such an effort. 

At slightly higher temperatures than 1 K, the available 
evidence is also due to [99] . The data concern a mixture 
of superfluid and normal helium and it is hard to dis- 
entangle the two. The measurements of [101], though 
intrinsically exciting in addition to having instigated the 
recent recent interest in the problem, are indirect. Here, 
one measures the decay of superfluid vorticity (with its 



own caveats) and notes that the behavior is similar to 
that of the classical vorticity. From this one can com- 
pute the energy dissipation rate and infer the classical 
Kolmogorov spectrum. 

Our conclusion is not that the —5/3 power is ruled 
out, but that the evidence is soft at present; one needs 
to produce more direct and convincing evidence. 

There is another interesting wrinkle. If one assumes 
that the wavelength of the Kelvin waves which dissipate 
or radiate the energy are very small compared to the Kol- 
mogorov scale, it is plausible to infer the spectral ampli- 
tude of fluctuations of superfluid velocity in the inter- 
mediate region between Kolmogorov and Kelvin scales. 
Presumably, the only relevant parameter in that range 
is the strain rate at the Kolmogorov scale, quite like the 
situation of the passive scalar spectrum at high Schmidt 
numbers. It then follows from dimensional reasoning that 
one should expect a —1 power for the spectrum in that 
region. Although no one has measured the spectrum di- 
rectly, indications from the decay of superfluid vorticity 
are that the energy spectrum is consistent with a —3 
power law [106]. This behavior is poorly understood at 
present. 

B. Visualization of quantized vortex lines 

One of the exciting developments of recent few years 

is the visualization of quantized vortices and their re- 
connection using small neutral particles [94, 107]. These 
particles are made by the in situ freezing of mixtures of 
hydrogen and helium. While these visualization stud- 
ies have confirmed some interesting aspects of quantized 
vortices such as rings and reconnections. the particles are 
still too large compared with the diameter of the vortices 
(by a factor of about 10'*). Thus, while it is easy to con- 
vince oneself that the particles get attracted to vortex 
cores and decorate them, it is obvious that the parti- 
cles are not always passive. One can calculate conditions 
under which the inertia of the particles has marginal in- 
fluence on vortex lines, but there is no controlled means 
to ensure that this happens always: One would have to 
devise smaller particles before one can be confident of the 
fine details. 



C. Concluding remarks on superfluid turbulence 

At least in the initial stages when the study of super- 
fiuid turbulence was brought closer to classical turbu- 
lence community, one of the hopes was that one might 
be able to create enormous Reynolds numbers in modest- 
sized facilities using helium II. However, this has turned 
out to be impossible in principle (though, to avoid con- 
fusion, we should reiterate that this goal has been suc- 
cessfully reached with helium I and the gaseous phase). 
The difficulty with helium II is that the superfluid vor- 
ticity introduces an effective kinematic viscosity which is 
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of the same order as the kinematic viscosity of hehum I 
[101, 105]. While there is a lot to learn and understand 
about superfluid turbulence as a subject of intrinsic in- 
terest, it is unlikely that it will unlock new paradigm for 
understanding classical turbulence. 

The new directions of superfluid turbulence concern 
helium 3 at much colder temperatures. 

VIII. FINAL REMARKS 

If we are interested in discovering laws underlying sys- 
tems with many strongly interacting degrees of freedom 
and are far from equilibrium, it is important to begin 
with a study a few of them with the same rigor and con- 
trol for which particle physics, say, is well known. We can 
probably make the case that hydrodynamic turbulence, 
which arises in flowing fluids, is an ideal paradigm. Our 
first point is that the dynamical equations for the motion 
of fluids are known to great accuracy, which means that 
understanding their analytic structure can greatly sup- 
plement experimental queries; in just the same way, com- 



puter simulations — even if they require much investment 
of time and money — can be far more useful here than for 
many other problems of the condensed phase, in which 
the interaction potential among microscopic parts is of- 
ten simply an educated guess. The stochasticity of tur- 
bulence (and of all systems that are driven hard) means 
that one may discern only laws that concern statistical 
behavior. If we are fortunate, these laws are universal in 
some well-understood sense. This is the way we regard 
the "problem of turbulence" . 

While we have not yet reached a state when we can 
declare victory (perhaps that may never happen in a 
strict sense), the "problem of turbulence" is being slowly 
chipped away by understanding, albeit partially, its sev- 
eral aspects. This review has touched a few aspects of the 
problem in which considerable progress has been made 
recently. There is, of course, much to do, and one needs 
to understand both the richness of the problem and the 
discipline needed to make a dent in one of its non-trivial 
aspects. 
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